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Popular Summary 

One of the primary instruments on both the Earth Observing System (EOS) Terra and 
Aqua satellites is the Moderate Resolution Imaging Spectroradiometer (MODIS). Using 
thermal-infrared bands from MODIS and heritage algorithms, we have developed an 
algorithm that determines the surface temperature of sea ice, known as the ice-surface 
temperature (1ST) under clear skies at 1-km resolution on a daily, global basis. In this 
paper, we present the MODIS 1ST algorithm and show two case studies of validation for 
the algorithm. One case study, which consisted of 255 samples from MODIS scenes over 
the South Pole, compared air temperature from the South Pole station with MODIS- 
derived ISTs. Results showed an uncertainty of <1.7°K in the MODIS measurement. A 
second case study consisted of 25 samples from MODIS scenes that were compared with 
air temperatures derived primarily from drifting buoys in the Arctic Ocean, and these 
results showed an uncertainty of <1.6°K in the MODIS measurement. In both cases, the 
air temperatures were acquired at a height of about 2 m or greater. We show that the 
uncertainty in the MODIS measurements would be even lower if actual ISTs could be 
compared with the MODIS-derived ISTs (since the air temperature at a 2-m height above 
the ground or ice is different from the temperature of the actual surface). Sea ice 1ST is 
an important quantity that is necessary for energy-balance modeling. The MODIS 1ST 
maps are the highest-resolution global, daily 1ST maps available. 
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Abstract 

Global sea ice products are produced from the Earth Observing System (EOS) Moderate 
Resolution Imaging Spectroradiometer (MODIS) on board both the Terra and Aqua 
satellites. Daily sea ice extent and ice-surface temperature (1ST) products are available at 
1- and 4-km resolution. Validation activities have been undertaken to assess the accuracy 
of the MODIS 1ST product at the South Pole station in Antarctica and in the Arctic Ocean 
using near-surface air-temperature data from a meteorological station and drifting buoys. 
Results from the study areas show that under clear skies, the MODIS ISTs are very close 
to those of the near-surface air temperatures with a bias of -1 . 1 and -1 .2°K, and an 
uncertainty of 1.6 and 1.7°K, respectively. It is shown that the uncertainties would be 
reduced if the actual temperature of the ice surface were reported instead of the near- 
surface air temperature. It is not possible to get an accurate 1ST from MODIS in the 
presence of even very thin clouds or fog, however using both the Advanced Microwave 
Scanning Radiometer-EOS (AMSR-E) and the MODIS on the Aqua satellite, it may be 
possible to develop a relationship between MODIS-derived 1ST and ice temperature 
derived from the AMSR-E. Since the AMSR-E measurements are generally unaffected 
by cloud cover, they may be used to complement the MODIS 1ST measurements. 


Introduction 

The presence of sea ice influences the temperature and circulation patterns of both the 
atmosphere and the oceans. Sea ice reduces the amount of solar radiation absorbed at the 
ocean surface and, with the overlying snow cover, serves as an insulator, restricting 
exchanges of heat, momentum and chemical constituents between the atmosphere and the 
ocean. Its large area coverage, 5-8% of the ocean surface, makes the sea ice cover a key 
parameter in the Earth’s energy balance. 

The importance of obtaining long-term, climate-data records of sea ice in the form of 
high-resolution, validated global sea ice maps, is highlighted by recent studies showing 
decade-scale changes in the global sea ice extent. Several authors have shown that the 
extent of sea ice in the Arctic has decreased since 1979 by ~3% per decade [l]-[3] and in 
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addition, the perennial sea ice cover in the Arctic declined at ~9% per decade from 1978 
to 2000 [4]. In the Southern Ocean the extent of sea ice increased ~1% per decade from 
1979-1998 [5]. Most of these studies have used space-borne coarse-resolution (~25 km) 
passive-microwave instruments to monitor the sea ice. 

Two classes of sensors, microwave and multispectral radiometers, are typically used for 
global mapping of sea ice extent, concentration, type and temperature. The ability of 
microwave instruments, both passive and active, to collect data through cloud cover and 
polar darkness makes them well suited to global monitoring of sea ice extent and 
dynamics. Passive-microwave data can also provide a measure of ice temperature [6]-[9] 
though the location within the snow/ice that the temperature represents varies with ice 
type. However, microwave instruments cannot collect data on sea ice albedo, needed for 
analysis of energy exchange between the ocean and atmosphere, and they have a 
relatively coarse spatial resolution which limits their utility for detailed studies of sea ice 
dynamics. In contrast, visible and/or near-infrared sensors can provide detail on the ice 
extent and concentration during clear-sky conditions as well as snow/ice albedo and ice- 
surface temperature (1ST). The Moderate Resolution Imaging Spectroradiometer 
(MODIS), flown on both the Terra and Aqua satellites (launched in December 1999 and 
May 2002, respectively), is useful for determination of sea ice extent, 1ST, albedo, 
movement, type, and concentration at resolutions of 250 m - 1 km. 

Daily extent and 1ST maps are produced using MODIS data at 1-km resolution. In this 
paper, we describe the algorithm used to create the MODIS 1ST maps, and provide 
validation information for this product. 


Background 

Global year-round records of ice-covered Arctic and Antarctic seas were acquired 
from several satellite passive-microwave radiometers beginning with the Nimbus-5 
Electrically Scanning Multichannel Microwave Radiometer (1972-1977), followed by the 
Nimbus-7 Scanning Multichannel Microwave Radiometer (1978-1987), and then by the 
series of Defense Meteorological Satellite Program Special Sensor Microwave/ Imagers 
(1987-present). The acquisition of microwave radiances from these sensors has allowed 
the production of hemispheric sea ice concentration maps and global time series of sea 
ice extent and area [1], [7], [10] & [1 1]. Sea ice maps have also been available from the 
National Oceanic and Atmospheric Administration (NOAA) National Ice Center (NIC) 
since 1972 [12] and are currently available from MODIS (from Febmary 24, 2000, to the 
present). 

Measurement of sea-ice 1ST is possible with optical sensors such as the Advanced Very 
High Resolution Radiometer (AVHRR), Landsat Thematic Mapper (TM) and Enhanced 
Thematic Mapper Plus (ETM+) [13]-[16] and MODIS. Comparisons of AVHRR-derived 
ice surface temperature with measurements during the Surface Heat Budget of the Arctic 
(SHEBA) experiment [17] and other field experiments show that the clear-sky 1ST can be 
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estimated with an accuracy (root-mean square error) of 0.5 - 1.5°K, depending on the 
season [18], [19] & [15]. 

For the retrieval of clear-sky 1ST, a split-window technique is used, where “split- 
window” refers to brightness temperature differences in the 1 1-12 pm atmospheric 
window. This technique allows for the correction of atmospheric effects, primarily water 
vapor. First employed to determine sea-surface temperature (SST) [20], the technique 
was subsequently used to determine 1ST in the Arctic with the AVHRR on NOAA polar- 
orbiting satellites [14]. This work was later refined and expanded to include ice surfaces 
in both polar regions and a broader range of atmospheric conditions. Several other 
investigators have employed variations of the split-window method for estimation of 1ST 
in the polar regions using the AVHRR and other instruments (e.g., [13], [21] & [22]). 

MODIS is an imaging spectroradiometer that provides imagery of the Earth’s surface and 
clouds in 36 discrete, narrow spectral bands from approximately 0.4 to 14.0 pm [23]. 
Global vegetation and land cover, global land-surface change, vegetation properties, 
surface albedo, land and ice surface temperature and snow and ice cover [24] products 
from MODIS are available on a daily or regular basis [25]. The spatial resolution of the 
MODIS instrument varies with spectral band and ranges from 250 m to 1 km at nadir. 


MODIS sea ice algorithms 

Algorithms for deriving the sea ice map products from MODIS are similar for both the 
Terra and the Aqua MODIS instruments, and are used to generate maps of sea ice extent 
as determined by reflectance and 1ST [26] & [27]. The products are produced globally at 
both 1-km and ~4-km spatial resolution. During daylight hours, two types of sea ice map 
products are produced; one is derived from reflectances and the other from 1ST while at 
nighttime, only the 1ST is produced (Table 1). MODIS sea ice products are generated 
automatically using MODIS sensor radiance data [28], a geolocation product [29], and 
the MODIS cloud mask product [30], [25] & [26]. In addition to sea ice extent and 1ST, 
quality assurance, latitude and longitude, and other information is also included with the 
data products. 

Fig. 1 shows the swath coverage acquired in a region of the Arctic Ocean north of Alaska 
from Level IB (LIB) data (left image) and the 1ST map product (right image). LIB is a 
swath (scene or granule) of MODIS data geolocated to latitude and longitude centers of 
1-km resolution pixels and is available in each of the 36 bands. The MODIS LIB image 
(left) is shown as a composite of MODIS bands 1 (0.62-0.67 pm), 4(0.545-0.565 pm) and 
3 (0.459-0.479 pm). The MODIS 1ST map (right) also shows ice extent because any 
pixel that has a temperature of less than 271.5°K is considered to be ice, and any pixel 
with a temperature of 27 1 ,5°K or greater is considered to be open water. Because sea ice 
is saline, it freezes at a temperature that is less than 273° K; we use 271.5° K as the cutoff 
temperature between water and ice, although a user can select his/her own cutoff 
temperature based on the 1ST to develop a sea ice extent map using the 1ST. 
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MODIS sea ice data products are produced as a sequence of products, beginning with a 
swath (scene) at a nominal pixel spatial resolution of 1 km and a nominal swath coverage 
of 2330 km (cross track) by 2030 km (along track) (about five minutes of satellite travel), 
and ends with a global-daily map product (Table 1). Automated selection of the most 
favorable observation from all the swaths acquired during a day generates a daily gridded 
sea ice product. A daily global product of sea ice extent and 1ST at ~4-km resolution in 
the Equal- Area Scalable Earth Grid (EASEGrid) projection [3 1] 
http://nsidc.org/data/ease/ has been developed and is being tested (Fig. 2). Production of 
this product is expected in the fall of 2003. 

There are several different data-product levels starting with Level 2 (L2). An L2 product 
is a geophysical product that remains in latitude and longitude orientation; it has not been 
temporally or spatially manipulated. A level-2G (L2G) product is a gridded format of a 
map projection. At L2G the data products are referred to as tiles, each tile being a piece, 
e.g., 10° latitude by 10° longitude area, of a map projection. Level-2 data products are 
gridded into L2G tiles by mapping the L2 pixels into cells of a tile in the map projection 
grid. The L2G algorithm creates a gridded product necessary for development of the 
level-3 (L3) products. An L3 product is a geophysical product that has been temporally 
and or spatially manipulated, and, for sea ice, is gridded into the final product, also an L3 
product, a daily global sea ice product for the north and south polar areas in the EASE- 
Grid projection at ~4-km resolution (Fig. 2). 

Both the Terra and Aqua MODIS products may be ordered through the National Snow 
and Ice Data Center (NSIDC) Distributed Active Archive Center (DAAC) Earth 
Observing System (EOS) Data Gateway (EDG) at 

http://nsidc.org/~imswww/pub/imswelcome/index.html T321. MOD29 is the Earth 
Science Data Type (ESDT) of the MODIS sea ice products for Terra, and MYD29 is the 
ESDT for Aqua. “Collection 4” or Version 4 represents the most advanced algorithm, 
and the products in Collection 4 are improved relative to Collection 3 and earlier 
collections although at the time of this writing, Collection 4 reprocessing is ongoing but 
incomplete. When Collection 5 begins in early 2004, all the data will be reprocessed 
beginning February 24, 2000, to the present into a consistent data set. 

The content of the sea ice data products differs between day (MOD29P1D or 
MYD29P1D) and night (MOD29P1N or MYD29P1N) because MODIS visible data are 
not acquired when the sensor is in night mode. The day product (MOD29P1D or 
MYD29P1D) contains both day and nighttime sea ice extent by reflectance and extent by 
1ST while the night product (MOD29P1N or MYD29P1N) contains only 1ST data 
processed from MODIS thermal data acquired during the night mode operation of the 
sensor [26]. These designations should be used when ordering MODIS data from 
NSIDC. 

Sea ice surface temperature. The 1ST algorithm is the same for both the Terra and Aqua 
MODIS. For 1ST, the split-window technique is implemented as a simple regression 
model of the form 
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T s =a + bT n +c(T n -T n ) + d[(T u - T n ){ sec 0-1) [1] 

where T s is the surface temperature, T\\ and Tn are the satellite measured brightness 
temperatures in the 1 1 and 12 pm channels (MODIS channels 31 and 32), 0is the sensor 
scan angle, and a, b, c, and d are regression coefficients. The split-window temperature 
difference is proportional to the atmospheric water vapor amount, and the scan angle 
provides information on the path length. 

To determine the empirical relationship in eq. [1], radiosonde data from drifting ice and 
land-based stations were used with a radiative transfer model, LOWTRAN [33], to 
simulate the sensor brightness temperatures. The sensor spectral response functions were 
convolved with the calculated radiances and converted to brightness temperatures for 
each radiosonde profile. The surface temperatures used in the model calculations were 
then regressed against the simulated brightness temperatures to determine the coefficients 
in eq. [1]. Radiosonde data poleward of 65° latitude were taken from three archives: the 
North Pole archive, the National Center for Atmospheric Research (NCAR) archive, and 
the Historical Arctic Rawinsonde archive [34]. Profile data for Antarctica are from the 
Antarctic Radiosonde Data Set [35], a compilation of data from 18 land stations, 
including 16 coastal stations and two interior stations. More than 1000 profiles were 
used in the each of the Arctic and Antarctic analyses. The spectral, scan angle-dependent 
emissivities of snow as specified in [15] were used in the radiative transfer calculations. 
The regression coefficients in eq. [1] were determined separately for the Arctic and 
Antarctic and for three temperature ranges. 

For the AVHRR, some SST methods also employ the 3.7 pm channel. Of course, the use 
of this channel would be limited to nighttime analyses because it measures emitted 
thermal energy as well as reflected solar radiation. However, this channel is often noisy 
at low temperatures when the amount of emitted energy is small, typical of the polar 
night. For the AVHRR it was found that the inclusion of the 3.7 pm channel in eq. [1] 
reduced the regression root-mean-square error by only 0.003 K, which we do not 
consider a significant improvement in accuracy. Nevertheless, both the 3.7 pm and 3.9 
pm MODIS channels (channels 20 and 21) will be investigated further in the future to 
determine if their use will improve the MODIS-derived 1ST. 

Fig. 3a shows a MODIS 1ST map of sea ice in the northern Greenland Sea. Note the 
lower ISTs in this map as compared to those in the Arctic Ocean north of Alaska (Fig. 1). 
The climatological surface air temperature is colder by more thanl0°K north of 
Greenland as compared to the Arctic Ocean just north of Alaska (see data from R. 
Jenne/NCAR in [1 1], In both scenes, the recently-refrozen leads are obvious as 
evidenced by their higher relative temperatures. Well-developed ice floes are evident in 
the Greenland Sea (Fig. 3a) because this is an active area and near the ice edge, while the 
area north of Alaska (Fig. 1) is a much more stable area and thus the ice has not been 
broken apart into distinct floes. 

Note the area of the highest ISTs in Fig. 3a, shown in the various shades of red in the 
southeastern part of the scene. This corresponds to an area of “warmer ice” or of very 


6 


low ice concentration approaching the ice edge. Open water (with an 1ST exceeding 
271.5°K) is also visible there. Though the ice edge cannot be seen on the full swath of 
the MODIS 1ST product due to obscuration by cloudcover, it is evident on an 89 GHz 
vertically-polarized Advanced Microwave Scanning Radiometer-EOS (AMSR-E) image 
shown in Fig. 3b. The location of the MODIS image (Fig. 3a) is shown in the box east- 
northeast of Greenland in Fig. 3b. The scale represents microwave brightness 
temperatures (T B s) from 180 - 300°K. Note that there is a general trend toward lower T B s 
(Fig. 3b) and lower ISTs (Fig. 3a) toward the northwest in the rectangular study area. 

Land and Cloud Masking in the Data Products. Land is masked using the 1 -km 
resolution United States Geological Survey (USGS) global land/water mask [36] stored 
in the MODIS geolocation product [29]. The land/water mask is used by both the sea ice 
algorithm and cloud-masking algorithm to determine if a pixel is land or ocean for 
processing in the algorithms. An issue with this land/water mask is that ice shelves, such 
as the Ross Ice Shelf of Antarctica, are classified as land and not as part of the ocean 
because they emanate from the ice sheet. Consequently, ice shelves have not been 
included in the MODIS sea ice products. Improved coastline data that maps ice shelves 
as ocean have been developed and will be integrated into a future improved land/water 
mask. Following its implementation, it will be possible to monitor large icebergs that 
calve from the ice shelves in Antarctica. 

Many situations occur that make it very difficult to identify clouds over sea ice in the 
daytime or nighttime even with the many visible and infrared MODIS bands. If the 
MODIS cloud mask reports a “certain cloud” condition for a pixel, then the sea ice 
product is reported as cloud for that same pixel [30]. If the cloud mask reports that the 
observation is “probably clear,” or “confident clear,” then it is interpreted as clear and the 
sea ice algorithm analyzes the pixel for 1ST. 

The accuracy of the MODIS cloud mask over snow and ice varies, with the highest 
accuracy generally found for daytime data and the lowest accuracy in darkness. Low 
illumination and the similar reflectance characteristics of clouds and sea ice at visible 
wavelengths cause difficulties in cloud masking using MODIS data. Additionally, 
atmospheric inversions can cause surface temperatures to be colder than the temperature 
of some clouds, e.g., [37], Snow/cloud discrimination problems are compounded during 
the polar night because the MODIS reflective bands are not useful during darkness. The 
absolute accuracy of the MODIS cloud mask used in the sea ice algorithm is not known. 
The Aqua MODIS cloud mask algorithm was revised to use band 7 (2. 1 pm) in place of 
band 6 (1.6 pm ) because most of the detectors on the Aqua MODIS in band 6 are non- 
functional. The revised Aqua cloud mask algorithm has been in operation since May 20, 
2003, and has affected the masking of clouds over sea ice relative to the cloud mask used 
with the Terra MODIS products. 

A particular problem for mapping clouds over sea ice is the identification of thin clouds 
and fog. Fog often develops in sea ice covered waters when water vapor forms following 
the formation of leads and polynyas. While it may be possible to see the sea ice on a 
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MODIS image, and even map the sea ice through the fog, it is not possible to retrieve an 
accurate 1ST through fog, even with the many visible and infrared MODIS bands. 


1ST Validation 

In order to determine the accuracy of the MODIS-derived 1ST measurements, MODIS- 
derived skin temperature and the observed surface air (not skin) temperature at South 
Pole were compared (Fig. 4). Though a land site, the South Pole is a good study site 
because surface-based cloud measurements are available. The results represent primarily 
nighttime cases, and are from April through December 2001 . Scenes were identified as 
clear by examining micropulse lidar data at the South Pole station. The closest pixel was 
used, i.e., no averaging of adjacent pixels was done. There is a bias -1.2°K, and a root- 
mean-square error (RMSE) or uncertainty of 1 .7°K where the satellite-derived skin 
temperature is, on average, lower than the measured surface air temperature. This 
difference is almost certainly due to the skin and air temperature differences that are 
common in the presence of atmospheric temperature inversions. The variability in the 
two measurements is greatest for low surface temperatures, possibly due to the lower 
radiometric accuracy of the MODIS infrared channels at such extreme temperatures. 

Validation in the Arctic Ocean was also conducted. Data from the National Oceanic and 
Atmospheric Administration (NOAA) Center for Operational Oceanographic Products 
and Service (NOAA/NOS CO-OPS) program 

(http://www.co-ops.nos.noaa.gov/co-QDs.htmn meteorological station located at the 
Prudhoe Bay tide station data in Prudhoe Bay, AK (70° 24.0' N, 148° 3 1 .6' W) were 
studied. The Prudhoe Bay tide station, established in July 1990, is located at the end of a 
causeway north of Prudhoe Bay. The air temperature measurements are acquired at a 
height of ~9 m above sea level [38]. The tide station also provides information on tides 
(water level) and barometric pressure, wind speed, direction and gust, including hourly 
air and water temperature data. 

An average of MODIS pixels nearest to the NOAA tide station at Prudhoe Bay, Alaska, 
was used to determine the 1ST at the time of the MODIS overpass. The pixel nearest the 
Pmdhoe Bay tide station was used, along with 2 to 5 immediate surrounding pixels 
(exclusive of cloud and land pixels). This was matched with the nearest (in time) air 
temperature provided by the NOAA/NOS CO-OPS program through the website shown 
above. 

Air temperatures from drifting buoys were also used to validate the Terra and Aqua 
MODIS 1ST maps. The North Pole Environmental Observatory (NPEO) provides long- 
term measurements that are needed to track changes in the Arctic environment, including 
an automated drifting station consisting of clusters of buoys fixed to the drifting sea ice 
[39]. The buoys, tracked with global-positioning system receivers, measure 
atmospheric, ice and other geophysical quantities in the Arctic Ocean. Air temperature is 
acquired at a height of ~1 .2 m. The Japan Marine Science and Technology Center 
(JAMSTEC) Compact Arctic Drifter #1, or J-CAD 1, was the first in a series that began 
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acquiring measurements in April 2000 in the Arctic Ocean. For validation of the MODIS 
1ST we used data from J-CAD 4 and 5 as well as from the Prudhoe Bay tide station. 

Radionov et al. [40] show that in the Arctic Basin during the winter when the cloud cover 
is continuous, the snow surface temperature did not differ from the air temperature 
measured at a height of 2 m. However, under cloud- free conditions, the difference 
between those temperatures reached 6-7°K, the snow temperature being sharply lower 
compared with the 2-m air temperature due to longwave radiative cooling when the sky 
was clear. During the Arctic cold period, defined as October through May, the average 
temperature of the snow (250 measurements per month) was lower by an average of 
0.6°K than the 2-m air temperature. During the warm period, defined as June through 
August, the average temperature of the snow was 0.7°K lower than the 2-m air 
temperature (125 - 246 measurements per month) [40]. Based on this information, a bias 
of at least 0.6°K should be considered when comparing air temperature measurements 
from the Prudhoe Bay Co-op station and the drifting buoys to the MODIS 1ST retrievals. 
All of our data (Table 2) were acquired during the “cold period” as defined above. Since 
the Prudhoe Bay in-situ air temperatures were acquired from a height of ~9 m, one would 
expect this difference to be even greater. 

Table 2 summarizes the Prudhoe Bay tide station and the Arctic drifting buoy results. 

Fig. 5 shows a comparison of the meteorological station or buoy temperatures and the 
MODIS-derived ISTs using 25 scenes that were deemed cloud-free at the measurement 
site. This data set reveals a bias of-l.l°K and aRMSE or uncertainty of 1.6°K, values 
that are close to those calculated for the South Pole data set (Fig. 4). If we adjust the air 
temperatures of the co-op station and the drifting buoys in accordance with the Radinov 
et al. results for the “cold period,” 0.6°K, then the adjusted bias is -0.46°K with an 
uncertainty of 1.3°K. However, because the meteorological station data and buoy data 
were acquired at a height of 1.2 and 9 m, respectively, and because the 1ST retrievals are 
for clear-sky conditions, it is difficult to have confidence in that exact number since the 
Radinov et al. measurements were made at a height of 2 m. 


Limitations 

A considerable amount of visual study was required to avoid selecting MODIS pixels that 
appeared to be cloud contaminated. Only about 14% of the MODIS scenes selected were 
used. The rest of the scenes were rejected for a variety of reasons. The main reason for 
rejection of a daytime scene was that clouds were observed in the LIB companion scene 
near the buoy or tide station. Additionally, when the difference between the MODIS and 
in-situ 1ST was large (~3° or greater) and there was a possibility of cloud contamination 
(e.g., clouds were nearby), or there was no LIB scene available due to darkness, the 
scene was rejected. As a result, the results in Table 2 are considerably better than would 
be the case if the scenes were selected randomly from “cloud-free” parts of the MODIS 
maps as determined from the MODIS cloud mask. 
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Though the MODIS-derived 1ST algorithm performs extremely well under totally clear- 
sky conditions, polar weather is often a factor that influences negatively the quality of the 
derived 1ST. It is virtually impossible to determine whether conditions are completely 
cloud free when deriving 1ST in an automated way using the MODIS cloud mask 
algorithm summary result flag of “certain cloud.” That summary result is based on 
several cloud tests applied in the cloud mask algorithm and confidences associated with 
the tests. In some situations the tests can fail to detect clouds or may detect clouds falsely 
resulting in an inaccurate result. Improved cloud masking results have been attained in 
the MODIS snow algorithm when selected cloud tests from the MODIS cloud mask are 
used to create a cloud mask specific to snow-covered land surfaces. Selective use of 
individual cloud tests from the MODIS cloud mask product should result in improved 
cloud masking relevant to the 1ST algorithm as has been the case for the MODIS snow 
maps [41]. 


Discussion and Conclusions 

The 1-km MODIS daily 1ST product represents the only global 1ST product available at 
such a fine resolution and may be used for energy-balance modeling. The clear-sky 
limitation can possibly be obviated using passive-microwave data in conjunction with 
MODIS data. Passive-microwave-derived ice temperatures are not directly comparable 
with ISTs because the passive-microwave temperature is an integrated measure of the 
temperature of the upper layers of the sea ice or/and overlying snow cover and the depth 
of penetration depends on the emissivity of the snow/ice. However, if a relationship can 
be determined between MODIS-derived ISTs and passive-microwave-derived ice 
temperature, then the passive-microwave derived ice temperatures can be used to 
approximate 1ST even through clouds. Such a product would provide modelers global 
1ST data daily through cloudcover and darkness though at a coarser resolution than is 
possible from MODIS if MODIS could be used alone. 

With both the AMSR-E and the MODIS on the Aqua satellite, it will be possible to assess 
the relationship between MODIS-derived 1ST and ice temperature derived from the 
AMSR-E. The recent launches of the AMSR sensors on the EOS Aqua and Advanced 
Earth Observing Satellite II (ADEOS II) spacecraft and the planned Defense 
Meteorological Satellite Program (DMSP) Special Sensor Microwave Imager/Sounder 
(SSMIS) and the National Polar-orbiting Operational Environmental Satellite System 
(NPOESS) missions will extend the sea ice records into the foreseeable fixture. 

The MODIS sea ice products represent a suite of products that is useful for long-term 
climate-data records. Additional work will need to be done in order to validate the 
products for a variety of ice types and geographic locations. However the MODIS 1ST 
product appears to be of comparable or superior quality as compared to predecessor data 
sets, such as products derived from the AVHRR, and should compare well with products 
that will be derived from future sensors for sea ice extent and 1ST. 
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The MODIS-derived 1ST provides an excellent measurement of the actual temperature of 
the surface of the sea ice under clear skies. Two different validation studies from the 
South Pole and the Arctic Ocean show excellent agreement, with a bias of - 1 . 1 and - 1 .2 
°K, and an uncertainty of 1.6 and 1.7 °K, respectively. 
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Figure 1. Left - True-color MODIS image using bands 1 (0.62-0.67 pm), 4(0.545-0.565 
pm) and 3 (0.459-0.479 pm) derived from Terra MODIS Level IB data, March 6, 2003 
(22:45 UTC), in the Arctic Ocean north of Alaska; Right - MODIS ice-surface 
temperature (1ST) map product of the same area at the same date and time. The location 
of the J-CAD 5 drifting buoy is shown inside the small yellow box at 76.91° N, 162.50° 
W (upper left). The color scale represents 1ST in °K. 

Figure 2. Prototype of the sea ice-surface temperature (1ST) global EASE-Grid 
composite product of the north polar area (average of May 15-19, 2000), MOD29C2. 
Color key shows ISTs from 230-27 1°K. The color scale represents 1ST in °K. 

Figure 3. a) Terra MODIS ice-surface temperature (1ST) product (MOD29) acquired on 
March 12, 2003 (18:45 UTC), in the northern Greenland Sea in the Arctic Ocean. The 
approximate center point of the image is 81.7° N, 1.0° E. The color scale represents 1ST 
in °K. b) Advanced Microwave Scanning Radiometer-EOS (AMSR-E) 89 GHz 
vertically-polarized image from the Aqua satellite, averaged from several swaths of data 
acquired on March 12, 2003, at a spatial resolution of ~5 km. The rectangular box east- 
northeast of Greenland is the area of the MODIS 1ST image from Fig. 3a. Microwave 
brightness temperatures (T B s) from 180-300°K are shown in the scale. 

Figure 4. MODIS 1ST retrieved temperatures and measured surface temperatures at the 
South Pole in Antarctica. Results are from a comparison of MODIS-derived skin 
temperature and the observed surface air (not skin) temperature at South Pole, primarily 
nighttime cases for the period April-December 2001. Scenes were identified as clear by 
examining micropulse lidar data at South Pole. The closest pixel was used, i.e., no 
averaging was done. 

Figure 5. MODIS 1ST retrieved temperatures and measured surface temperatures at three 
locations in the Arctic Ocean on various dates between January 25, 2002, and May 9, 
2003. Most of the measured temperatures were acquired from buoy data (see text for 
explanation), but a few were acquired from the tide station at Prudhoe Bay, AK (see text 
for explanation). 
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Table 1. MODIS sea ice data products. 


Long Name 


Spatial Resolution 




MODIS/Terra Sea Ice Extent 5-Min 
L2 Swath 1km 

MOD29 

l-km resolution, swath of 
MODIS data 

MODIS/Terra Sea Ice Extent Daily L3 
Global EASE-Grid Day 

MOD29P1D 


MODIS/Terra Sea Ice Extent Daily L3 
Global EASE-Grid Night 

MOD29P1N 

l-km resolution, projected, 
gridded tile data 

MODIS/Terra Sea Ice Extent and Ice 
Surface Temperature Daily L3 Global 
4km EASE-Grid Day * 

MOD29E1D 

4-km resolution, global, gridded 

MODIS/Terra Ice Surface 
Temperature Daily L3 Global 4km 
EASE-Grid Night * 

MOD29E1N 

4-km resolution, global, gridded 

MODIS/Terra Sea Ice Extent and Ice 
Surface Temperature 8-Day L3 Global 
EASE-Grid Day ** 

MOD29E2D 

4-km resolution, global, gridded 

MODIS/Terra Ice Surface 
Temperature 8-Day L3 Global EASE- 
Grid Night** 

MOD29E2N 

4-km resolution, global, gridded 

MODIS/Terra Sea Ice Extent and Ice 
Surface Temperature Monthly 
L3GIobal 4km EASE-Grid Day** 

MOD29EMD 

4-km resolution, global, gridded 

MODIS/Terra Ice Surface 
Temperature Monthly L3Global 4km 
EASE-Grid Night ** 

MOD29EMN 

4-km resolution, global, gridded 

MODIS/ Aqua Sea Ice Extent 5-Min 
L2 Swath 1km 

MYD29 

l-km resolution, swath of 
MODIS data 

MODIS/Aqua Sea Ice Extent Daily L3 
Global EASE-Grid Day** 

MYD29P1D 

l-km resolution, projected, 
gridded tile data 

MODIS/Aqua Sea Ice Extent Daily L3 
Global EASE-Grid Night** 

MYD29P1N 

— ■ 

MODIS/Aqua Sea Ice Extent and Ice 
Surface Temperature Daily L3 Global 
4km EASE-Grid Day* 

MYD29E1D 

4-km resolution, global, gridded 

MODIS/Aqua Ice Surface 
Temperature Daily L3 Global 4km 
EASE-Grid Night * 

MYD29E1N 

4-km resolution, global, gridded 

MODIS/Aqua Sea Ice Extent and Ice 
Surface Temperature 8-Day L3 Global 
EASE-Grid Day ** 

MYD29E2D 

4-km resolution, global, gridded 

MODIS/Aqua Ice Surface 
Temperature 8-Day L3 Global EASE- 
Grid Night** 

MYD29E2N 

4-km resolution, global, gridded 

MODIS/Aqua Sea Ice Monthly Global 
EASE-Grid ** 

MYD29EMD 

4-km resolution, global, gridded 





♦Future enhancement or product (to be implemented in fall 2003). 


♦♦Future enhancement or product (to be implemented in 2004). 
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Table 2. 1ST validation summary for J-CAD 4 and 5 drifting buoys, and the 
Prudhoe Bay tide station. “Terra” and “Aqua” refer to images acquired from 
either the Terra or Aqua MODIS instruments, and 003 and 004 refer to either 
version 3 or 4 of the MODIS 1ST algorithm. Buoy or station temperatures are 
subtracted from the MODIS ISTs to derive the difference. 


J-CAD 4 - JAMSTEC buoy 


Date 

Time (UTC) - 
data specification 

MODIS 1ST 
°K 

J-CAD 4 buoy 
°K 

Difference 

°K 

Mar 1 2003 

1905 - Terra 004 

234.5 

235.9 

-1.4 

Mar 2 2003 

1810 -Terra 004 

239.1 

240.2 

-1.1 

Mar 3 2003 

0250 - Aqua 003 

237.1 

238.9 

-1.8 

Mar 3 2003 

0745 -Aqua 003 

237.8 

239.7 

-1.9 

Mar 3 2003 


238.3 

240.2 

-1.9 

Mar 9 2003 

1500 -Terra 004 

243.8 

243.9 



2105 -Terra 004 

269.7 

268.4 

1.3 


2010 -Terra 004 





J-CAD 5 - J AMSTEC buoy 


Date 

Time (UTC) - 
data specification 

MODIS 1ST 
°K 

J-CAD 5 buoy 
°K 

Difference 

°K 

Mar 6 2003 

2245 - Terra 004 

254.4 

256.7 

-2.3 

Mar 7 2003 

0200 - Terra 004 

252.1 

254.4 

-2.3 

Mar 29 2003 

0755 -Terra 004 

247.4 

249.9 

-2.5 

IJfflfJtWtkM 

0250 - Terra 004 

251.7 

253.9 

-2.2 


2240 - Terra 004 

250.6 

250.9 


ISSHE&SH 

2340 - Terra 004 

262.4 

262.4 



0255 - Terra 004 

261.5 

262.4 

-0.9 


Prudhoe Bay, AK - NOAA tide station 


Date 

Time (UTC) - 
data specification 

MODIS 1ST 
°K 

Prudhoe Bay °K 

Difference 

°K 

Jan 25 2002 

2230 - Terra 003 

235.3 

235.3 


Mar 3 2002 

2 100 -Terra 003 

245.0 

244.7 

IHE9H 

Mar 29 2002 

2145 -Terra 003 

248.1 

248.3 


Mar 30 2002 

2230 - Terra 003 

250.1 

248.5 

1.6 

Mar 31 2002 

2135 -Terra 003 

251.5 

251.1 



2205 - Terra 003 

272.9 

275.7 

-2.8 


21 10 -Terra 003 

273.3 

275.9 

-2.6 

Nov 20 2002 

2210 -Terra 003 

249.9 

252.5 

-2.6 

Jan 19 2003 

2055 - Terra 004 

254.9 

256.3 

-1.4 

Jan 19 2003 

2235 - Terra 004 

255.8 

257.4 

-1.6 













































































